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Abstract 

Calculations in unquenched QCD for the scalar glueball Spectrum have confirmed previous 
results of Gluodynamics finding a glueball at ~ 1750 MeV. I analyze the implications of this 
discovery from the point of view of glueball - meson mixing at the light of the experimental scalar 
sprectrum. 


Pacs: 12.38.-t, 12.38.Aw,14.70.-e, 14.80.-j 
Keywords: meson, glueball, mixing. 


1 Introduction 


Glueballs have not been an easy subject to study due to the lack of phenomenological support 
and therefore much debate has been associated with their properties [T]. I center the discussion 
here at the consequences of the Spectrum obtained by lattice QCD where the results seem to be 
converging. In quenched QCD the masses of the scalar glueballs appear large > 1700 MeV [2H5], a 
result which has been confirmed by unquenched calculations [6]. . 

Several glueball-meson mixing scenarios have been discussed in the literature using either lattice 
calculations or phenomenology [2l[3lE^9]. I implement a combination of lattice results and phe- 
nomenology to study the implications of recent lattice results in the possible mixing scenarios in 
the scalar sector. 

2 The scalar Spectrum 

Lattice QCD provides us with a value for the mass of the O"*"*" glueball states as shown in Table [TJ 



Mass MeV 

Unquenched 

Quenched 

Cl 

Mp 

Ky 

0++ 

1795(60) 

1730(50)(80) 

1710(50)(80) 

0++ 

3760(240) 

2670(180)(130) 



Table 1; Clueball masses with assignments. The column Cl reports the results of the un¬ 
quenched QCD calculation by Cregory et al. [6], the columns Mp and Ky show the data from 
Morningstar and Peardon [3j and Chen et al. [5] respectively. 

The three calculations give a similar mass for their lowest state for which I take the mean 
1743 ± 42 MeV in my analysis. 

In Table m I show the experimental scalar Spectrum, namely that of the particles labelled /q. 



Name 

Mass MeV 

Width MeV 

Comment 

0++ 

/o (500) 

400-550 

400-700 


0++ 

/o (980) 

990±20 

40-100 


0++ 

/o (1370) 

1200-1500 

200-500 


0++ 

/o (1500) 

1505 ± 6 

109 ± 7 


0++ 

/o(1710) 

1720 ± 6 

135 ± 8 


0++ 

/o (2020) 

1992 ± 16 

442 ± 60 

needs conhrmation 

0++ 

/o (2100) 

2013 ± 8 

209 ± 19 

needs conhrmation 

0++ 

/o (2200) 

2189 ± 13 

238 ±50 

needs conhrmation 


Table 2: Scalar particles appearing in the PDC’s particle listings m- 
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The three heaviest states have not been confirmed and some authors also question the existence 
of the /o(1370). 

An observation at the light of the experimental Spectrum is that the excited glueballs obtained 
in QCD calculations are very high in mass and therefore I do not expect them to mix with the 
mesons. Thus the lattice calculations and the observed scalar Spectrum lead to a scenario of one 
glueball amidst several scalar mesons in the range between 1 — 2 GeV. 


3 Glueball-Meson Mixing 

The most naive way of implementing a mixing scenario is by assuming a hamiltonian which is not 
diagonal in the unmixed states. Thus there are two ingredients to this hamiltonian, the unmixed 
masses and the off-diagonal mixing parameters. 

The calculations thus far have led to different realizations of mixing. Some obtain that the 
/o(1710) is mostly glueball with small admixtures of qq states [HEIIQ]. Others claim that the 
/o(1500) is mostly glueball, with small admixtures of qq states [HIH]. 



Figure 1: The largest eigenvalue of the studied analog of a three state mixing matrix is shown 
as a function of the relevant parameters in Cardano’s formulas. The valnes used for the diagonal 
elements were 1750,1500,1300 close to the physical /o masses, and the mixing parameters 6, c range 
from 0 —>• 1000, and a from —1000 —> 1000 in discrete steps. 


I analyze next which scenarios could arise in a full QCD calculation. The idea is to construct 
a hamiltonian which is not diagonal in the unmixed lattice states and diagonalize it to obtain the 
physical states. The decisive input is the glueball mass as given by lattice QCD. I assume that the 
rest of the particles in the unmixed Fock space are scalar mesons. Since the scalar meson Spectrum 
of lattice QCD is not well known I will be guided by phenomenology. 

For the off-diagonal matrix elements I shall use perturbation theory and 1/A^c arguments. With 
these assumptions the matrix elements are of the form ~ | < ^n\Hnon-diagonai\^m > \^/{En — Em), 
therefore they depend on the inverse of the decay coupling f meson ~ , fgiuebaii ~ -^c i on 

the inverse of the energy difference of the two levels. For almost degenerate unmixed states, the 
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mixing parameters might be larger than for states of quite different masses. Thus the hamiltonian 
matrix might look for a three state mixing as, 

^ nisi a 6 \ 

a ms 2 c , (1) 

\ b c mg ) 

where, and ms 2 represent the unmixed scalar meson masses, and mg the unmixed glueball mass. 
The diagonal matrix elements I take motivated by phenomenology and their values incorporate all 
higher order corrections, ~ 1 + 0{N^ ^) + ... for the glueballs, and ~ 1 + 0{N^ ^) + ... for the 
mesons. The mixing matrix elements are ~ 0{N~‘^) for meson-glueball states since Ai? ~ 
but they are also ~ 0{N~‘^) for the meson-meson states since A.E ~ 0(1) for them. Thus a,b and 
c will be of the same order of magnitude, i.e. < 250 MeV. 

It is trivial to show that in a two by two symmetric mass mixing matrix the the unmixed mass 
of the heavy state will increase after non trivial diagonalisation. For a three state mass mixing 
matrix it is not so trivial to show that the heaviest state mass also increases after non trivial 
diagonalisation. This can be done using Cardano’s formulas, as shown in Fig. [H Cardano’s 
formulas depend only of one sign, sign{abc), therefore one reaches all the desired values by only 
letting one of the parameters become negative. I have used as diagonal elements numbers close to 
the actual /o masses and the off diagonal elements have taken values far beyond my expectation 
for the physical mixing parameter values. 




Figure 2: Left: Eigenvalues of the mixing matrix for the glueball and two meson states. The masses 
and their er rors are from PDG m- Right; The probability distribution for the /o(1710) of the 
unmixed states as a function of c, having fixed a = b = 50 MeV. The dotted line represents the 
glueball with unmixed mass (1670 MeV), the full line the heavy meson with unmixed mass (1530 
MeV), the dashed line the light meson with unmixed mass (~ 1380 MeV). The error bars in the 
probability curves are just to guide the eye. 


The QCD lattice mass value for the scalar glueball is very high and the error relatively small if 
the three values shown in TableHJare averaged with errors in quadrature mg = 1743±42 MeV. One 
can still argue, given the errors in the lattice calculation, that the glueball mass could be below, 
but certainly not much below, the experimental /o(1710). Therefore, I foresee two scenarios: i) 
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the unmixed QCD mass value is below the /o(1710) but close to it, and ii) the unmixed glueball 
mass value is above the /o(1710) in agreement with the lattice QCD average value. In the latter 
case I will analyze two cases depending on the unmixed Spectrum, one will be associated with weak 
mixing, while the other with strong mixing. 




Figure 3: Left: Eigenvalues of the mixing matrix for the glueball and two meson states. The masses 
and their er rors are from PDG m- Right; The probability distribution for the /o(1710) of the 
unmixed states as a function of c . The other mixing parameters have been fixed to a = 100 MeV 
and b = 100. The dotted line represents the unmixed glueball (1750), the dashed line represents 
the light unmixed meson (1550), the full line represents the heavy unmixed meson (1900) 


The first scenario requires very small mixing, as can be seen in Figl2] (left). The unmixed valnes 
were chosen = 1380 MeV, ms 2 = 1530 MeV and mg = 1670 MeV, the latter two standard 
deviations below the central lattice value. The mixing parameters that fit the data are small 
a = 6 = 50 MeV and c = 70 MeV. Note that in this case one can also trivially construct a two state 
weak mixing scenario, with the outcome that the /o(1710) is again mostly glueball, the /o(1500) 
mostly meson and the mixing parameters are small.Thus the /o(1370) is not required from the 
point of view of mixing. 

I show in Figl2] (right) the probability distribution for the final mostly glueball state. It turns out 
that the experimental /o(1710) comes out mostly glueball in agreement with the lattice calculation 
of mixing mis] and other phenomenological analyisis [9]. 

Let me now elaborate on the other scenario, namely assuming mg > 1710 MeV, i.e. I take 1750 
MeV, close to the average value, as the unmixed glueball mass. As a result of Cardano’s analysis, I 
know that the only way to get that mass down to the /o(1710) is by incorporating a heavier meson. 
In particular, I show in Fig. [3] (left) a three state mixing where the /o(1710) is mostly glueball but 
the /o (2020), a not yet confirmed particle, must exist to get the experimental valnes for the lower 
masses. The corresponding probability curves, Fig. [3] (right), show that the glueball component is 
at the level of 70%. Most of the other 30 % is carried by the /o (2020). Resides the existence of 
the /o(2020), another feature of this scheme is that the /o(1370), whatever it be, is not required, 
and thus the error bars are really very small and the fit quite restrictive. Another characteristic of 
this weak mixing scheme in the spread distribution of the unmixed Spectrum; the particles appear 
with an energy step of 150 — 200 MeV. Note that in this case a two state mixing scenario could 
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Figure 4: Eigenvalues of the mixing matrix for the glueball and two meson states. The masses and 
their errors are from PDG m- 



b (MeV) 


also be trivially constructed, but here the /o(2020) would be unavoidable and the /o(1500) would 
decouple from the mixing scheme. 

Within the latter scenario also a strong mixing relization arises. As shown in FigØl I am able to 
get a good fit to the data starting from almost degenerate unmixed states. In this case the unmixed 
values are for two mesons at 1710 MeV and 1760 MeV and a glueball at 1750 MeV. The mixing 
parameters between the almost degenerate light meson and glueball is large b = 210 MeV, while 
the others are normal a = 100 MeV, c = 70 MeV. On the other hånd this strong mixing affects the 
/o(1710) strongly, which now contains very little glueball, ~ 17%, while increasing dramatically 
the glueball content of the /o(1500), ~ 41%, and some that of the /o(2020), ~ 42% , see FigO In 
this strong mixing scenario the glueball is a relic of the original glueball state and therefore it will 
be difficult to single out glueball properties. Again no need for the /o(1370) in the fit. 



b (MeV) b (MeV) b (MeV) 


Figure 5: The probability distributions as a function of b of the mixed states, for mixing parameters 
a = 100 MeV, c = 70 MeV. The notation follows the unmixed states: the full line represents the 
glueball (1750 MeV), the dashed line represents the meson (1710 MeV), the dotted line represents 
the meson (1760) MeV. The left figure represents the probability distribution of the /o(2020), the 
middle that of the /o(1710), and the right figure that of the /o(1500). 
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I might summarize at this stage my tindings by stating that accepting a mass valne for a 
glueball as obtained from lattice QCD for the unmixed glueball state and allowing for mixing 
under the strict scrutiny of the l/N^ expansion, two scenarios appear as dictated by the shape of 
the unmixed Spectrum: i) a week mixing scenario, whose consequence is that the /o(1710) is an 
almost pure glueball state; ii) a strong mixing scenario in which the /o(1710) has almost no glueball 
component, while the /o(1500) and the /o(2020) have about 40% glueball component. The role 
played by the scalar meson Spectrum in the realization of mixing is fundamental. A detailed lattice 
QCD description of this Spectrum would determine the mixing scenario. Within the weak mixing 
scenario it is very important the value of the mass of the unmixed glueball state. If it is below the 
/o(1710) the mixing proceeds via the /o(1500) and the /o(1370) and it is not very constrained due 
to the large indeterminacy in the latter. If that mass is above 1710 MeV the /o(2020) is required 
to bring the mass value down, while the /o(1370) is unnecessary. 

4 Conclusions 

In this paper I take a very pragmatic point of view. I assume that the lowest found lattice QCD 
glueball State is associated with real physical states within the /o Spectrum and try to establish the 
connection by implementing phenomenologically mixing with the nearby scalar mesons. In order 
to estimate the mixings I use the splitting pattern of the /o Spectrum and 1/A^c arguments. 

I find two possible weak mixing scenarios which associate the lattice glueball mostly with the 
/o(1710). In this case the exact value of the unmixed glueball state is fundamental to determine 
the physical realization of mixing, being the /o(1710) the dominant scale. If the unmixed glueball 
mass is below that value, the /o (1370) should exist and may play a relevant role. If the mass is 
above, at least the /o (2020) is required to enter the mixing to bring down the unmixed mass to the 
physical mass. The scalar Spectrum in both cases is loose, with an energy step between 150 — 200 
MeV. 

I also find a strong mixing scenario which leads to the /o(1500) and the /o(2020) having a strong 
glueball component, while the /o(1710) is mostly mesonic in character. In this case the unmixed 
scalar Spectrum presents an almost degeneracy around the unmixed glueball mass. 

The main conclusion of this calculation is that given the faet that mixings are difficult to calculate 
in lattice QCD, a good knowledge of the scalar meson Spectrum together with phenomenology will 
clarify the glueball constituency of the physical /o’s in full QCD. Once the scalar Spectrum is known 
the /o Spectrum and the arguments about the mixing used in this communication will fix quite 
strongly the mixing parameters. 

Giving the above analysis I find that the /o(1710) being mostly glueball is the most natural 
scenario. If it is accompanied, besides the /o(1500), by the /o(1370) or the /o(2020) in the mixing 
scenario is a matter of experimental determination. It is clear that for this purpose not only 
masses, but also decays will be needed to determine the main properties of the glueball component. 
In here I have presented a guide of possible mixing schemes as characterized by the structure of 
the Spectrum. 
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